THE N-ORDER ROGUE WAVES OF FOKAS-LENELLS EQUATION 



SHUWEI XU \ JINGSONG HE ^ *, YI CHENG ^ AND K. PORSEIZAN ^ 



^School of Mathematical Sciences, USTC, Hefei, Anhui 230026, P. R. China 
^Department of Mathematics, Ningbo University, Ningbo, Zhejiang 315211, P. R. China 
^Department of Physics, Pondicherry University, Puducherry 605014, India 

Abstract. Considering certain terms of the next asymptotic order beyond the nonhnear Schrodinger 
equation (NLS) equation, the Fokas-Lenells (FL) equation governed by the FL system arise as a model 
for nonlinear pulse propagation in optical fibers. The expressions of the g'"' and r'"' in the FL system 
are generated by n-fold Darboux transformation (DT), each element of the matrix is a 22 matrix, 
expressed by a ratio of (2n + 1) x (2n + 1) determinant and 2n x 2n determinant of eigenfunctions. 
Further, a Taylor series expansion about the n-order breather solutions q'"' generated using DT by 
assuming periodic seed solutions under reduction can generate the n-order rogue waves of the FL 
equation explicitly with 2n + 3 free parameters. 
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1. INTRODUCTION 

The Fokas-Lenells equation (FL) [IHS] 

iqxt - iqxx + 2gx - Qxqq* + iq = 0, (1) 

is one of the important models from both mathematical and physical considerations. In eq.([T|), q 
represents the complex field envelope and asterisk denotes its complex conjugation, the subscript x 
(or t) denotes partial derivative with respect to x (or t). The FL equation fP,^ is related to the 
nonlinear Schrodinger (NLS) equation in the same way as the Camassa-Holm equation associated 
with the KdV equation. The authors of [3j , after deriving associated Lax pair and using initial value 
problem, they were able to solve the equation. The soliton solutions of the FL equation have been 
constructed via the Riemann-Hilbert method in [2] and through dressing method in [5]. The breather 
solutions of the FL equation have also been constructed via a dressing-Backlund transformation related 
to the Riemann-Hilbert problem formulation of the inverse scattering theory [6]; also in the same 
paper three instability regions were analyzed, associated with a single unstable wave number. The 
FL equation actually describes the first negative flow of the integrable hierarchy associated with the 
derivative nonlinear Schrodinger (DNLS) equation j4] (it also belongs to the deformed DNLS hierarchy 
proposed by A.Kundu [ZIIH]). The lattice representation and the dark solitons of the FL equation have 
been presented in [9], where a relationship is also established between the FL equation and other 
integrable models including the NLS equation, the Merola-Ragnisco-Tu equations and the Ablowitz- 
Ladik equation. Recently, it has been shown that the periodic initial value problem for the FL 
equation is well-posed in a Sobolev space with exponent greater than 3/2 [10]. In optics, considering 
suitable higher order linear and nonlinear optical effects, the FL equation has been derived as a model 
to describe the femtosecond pulse propagation through single mode optical silica fiber and several 
interesting solutions have been constructed [4]. 
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Rogue waves have recently been the subject of intensive investigations in oceanography |llfll3] . 
where they occur due to either modulation instability [14H16] . or other processes |17fll9) . The first 
order rogue wave usually takes the form of a single peak hump with two caves in a plane with a 
nonzero boundary. One of the possible generating mechanisms for rogue waves is through the creation 
of breathers which can be formed as a consequence of modulation instability. Then, larger rogue waves 
can emerge when two or more breathers collide with each other |20fl22j . Rogue waves have also been 
observed in space plasmas p3H27] . as well as in optics when propagating high power optical radia- 
tion through photonic crystal fibers |28H30) . Though the rogue waves have been reported in different 
branches of physics where the system dynamics is governed by single equation, to the best our knowl- 
edge, they have been observed and reported very little in the coupled systems. For example, Rogue 
waves of the coupled Schrodinger equations are construct in the literature [3TH33] . In experiment, 
the rogue waves in a multistable system [M] is revealed by experiments with an erbium-doped fiber 
laser driven by harmonic pump modulation. Considering higher order effects in the propagation of 
femtosecond pulses, rogue waves have been reported in the Hirota equation |35p36j and in resonant 
erbium-doped fibre system governed by a coupled system of the nonlinear Schrodinger equation and 
the Maxwell-Bloch (NLS-MB) equations [37J. Very recently, the new types of matter rogue waves [38j 
have been reported in F = 1 spinor Bose-Einstein condensate governed by a three-component NLS 
equations. Some interesting results on the multi-rogue wave solutions of the NLS equation have also 
been done in |39H45j . which shows that there exist many patterns of the rogue waves and their for- 
mulae have extreme complexity. 

Considering the physical significance of the FL equation, inspired by the importance of the recent 
interesting developments in the rogue waves of the NLS type equations, so we have reported the first 
order rogue wave [46] of the FL equation in by Darboux transformation (DT) |47H49j . Our construc- 
tion reveals that there exists some deviations in their solutions and DT between the FL system and 
other integrable models such as Ablowitz-Kaup-Newell-Segur (AKNS) system [SOpSlj and Kaup-Newell 
(KN) system |52H54j . The purpose of this paper is to provide a detailed derivation of the determinant 
representation of the N-fold DT of the FL equation as we have done for the case of NLS equation [49] , 
which will be used to construct higher order rogue wave. Several different patterns of the higher order 
rogue waves will be plotted according to the determinant representation. 

The organization of this paper is as follows. In section 2, it provides a relatively simple approach 
to DT for the FL system, which is followed by the determinant representation of the n-fold DT and 
formulae of g'"' and rt""! expressed by eigenfunctions of spectral problem are given. The reduction 
of DT of the FL system to the FL equation is also discussed by choosing paired eigenvalues and 
eigenfunctions. In section 3, a Taylor series expansion about the n-order breather solutions generated 
by DT from a periodic seed solution with a constant amplitude can construct the n-order rogue waves 
of the FL equation in the determinant forms with 2n + 3 free parameters. Finally, we conclude the 
results in section 4. 

2. Darboux transformation 
Let us start from the non-trivial flow of the FL (Fokas-Lenells) system [3] , 

iqxt - kxx + - Qxqr + iq = 0, (2) 

irxt - irxx - '^r^ + r^rq + ir = 0, (3) 

which are exactly reduced to the FL eq.([T]) for r = q* while the choice r = —q* would lead to eq.([T]) 
with the sign of the nonlinear term changed. The Lax pairs corresponding to coupled FL eq.([2]) and 
([3]) can be given by the FL spectral problem [3l 

d^i; = (JA2 + QX)i; = Ui;, (4) 
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dt^ = (JA^ + Q\ + Vq + V-i\-^ + - JA-2)V' = V^, (5) 

with 

Here A is an arbitrary complex number, called the eigenvalue(or spectral parameter), and ij) is called 
the eigenfunction associated with A of the FL system. Equations ([2]) and ([3]) are equivalent to the 
integrability condition C/i - + [U, = of (g]) and ([5]). 

The main task of this section is to present a detailed derivation of the DT of the FL equation and 
the determinant representation of the n-fold transformation. Based on the DT for the NLS [47n49j 
and the DNLS |26 y 27 1[53ll54j . the main steps are : 1) finding 2x2 matrix T so that the FL spectral 
problem eq.dH and eq.([5|) is covariant, then get new solution (g[^],r[^l) expressed by elements of T 
and seed solution {q,r); 2) finding the expressions of elements of T in terms of eigenfunctions of FL 
spectral problem corresponding to the seed solution {q,r); 3) to get the determinant representation of 
n-fold DT Tn and new solutions (gl"], r^"!) by n-times iteration of the DT; 4) to consider the reduction 
condition: = (r^"!)* by choosing special eigenvalue A^ and its eigenfunction tpk, and then get g'"! of 
the FL equation expressed by its seed solution q and its associated eigenfunctions {^pk, k = 1,2, ■ ■ ■ ,n}. 
However, we shall use the kernel of n-fold DT(T„) to fix it in the third step instead of iteration. 

It is straightforward to see that the spectral problem (jj]) and ([5]) are transformed to 

= C/W , C/W = {T^ + T U)T-\ (6) 

= ^w, = (Tt + r v)t-\ (7) 

under a gauge transformation 

V^W =T tp. (8) 
By cross differentiating ([6]) and ([7]), we obtain 

f/W* - V^'K + [f/W, yW] = T{Ut -V, + [U, V])T-\ (9) 

This implies that, in order to make eqs.([2]) and ([3]) invariant under the transformation ([8]), it is crucial 
to search a matrix T so that [/[^l, have the same forms as U, V. At the same time the old 
potentiator seed solution) (g, r) in spectral matrixes U, V are mapped into new potentials(or new 
solution) (gt^l, A^^) in transformed spectral matrixes U^^\ V^^\ 

2.1 One- fold Darboux transformation of the FL system 

Considering the universality of DT, suppose that the trial Darboux matrix T in eq.([8]) is of the form 

where a_i, c_i, ao, bo, cq, do, ai,bi,ci, di are functions of x and t. From 

T^ + T U = [/W T, (11) 
comparing the coefficients of ,j = 3, 2, 1, 0, —1, it yields 
\^ :bi = 0, ci = 0, 

X'^ : qx ai + 2 i bo - q^^^^di = 0, -r^^K ai + di - 2 i co = 0, 
X^ : + rx bo- q^^^^co = 0, di^ + qa:Co - r^^^bo = 0, 
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Qxao - q^^^xdo + 2*6-1 = 0, -r^^^^ao + r^do - 2ic„i = 0, 
A° : -gWc_i + r^b-i + ao^ = 0, -g^d-i + q^a-i + 60^, = 0, 

- rWa_i + ra;(i„i + cq^ = 0, -rW6_i + q^C-i + = 0, 
A"^ : a-i,^ = b-i^ = c-i^ = d^i,^ = 0. (12) 

The terms in the previous equation a_i, c_i, are ah functions of t only. Similarly, from 

Tt + TV = T, (13) 

comparing the coefficients ofA-',j = 2,1,0,— 1,-2,— 3 under the condition 61 = ci = 0, it implies 

A"^ : 6-1 = c-i = 0, 

A"2 : + ga_i + 60 = 0, -rWa_i + r(i_i + cq = 0, 

A"^ : a_ii - -ig^co + -iq'WrWa-i - -iqra-i - -irbo = 0, -Q'Wdo + qao = 0, 

11 11 
d-it + -irWfto - -igWrWd_i + -igrd.i + -igco = 0, -rWao + rdo = 0, 

A° : - ^iqrao + ^igl^Vt^'ao = 0, dot + ^iqrdo - ^iq^'^^A^^o = 0, 

qxd-i + ^iqai - 2ibo - ^iq^^^di + 6ot - Q^^^xd-i + ^igW^Wfto + ^iqrbo = 0, 
rxd-i - ^irdi + 2ico + ^ir'^^ai + cot - r^^Ka-i - ^ig'^Vl^'co - ^iqrco = 0, 

A^ : -gW^do + g^ao = 0, -rl^^^ao + r^do = 0, 

^zgWrWai - ^iqrai + + r^bo - g^^co = 0, 

- ^igWr'^l^i + ^igrdi + dn + qxCo - r'^^^^o = 0, 
A^ : q^ai - ql^Ui + 2ibo = 0, r^di - g^oi - 2ico = 0. (14) 

In order to get the non-trivial solutions, we shall construct a basic (or one- fold) DT matrix T under 
an assumption that ao = and do = 0. If we set ao 7^ 0, then we can get that do is not zero by 
—q^^^do + qao = and — gt^'do + Qo-o = from eq. (fH|) . and furthermore find that some coefficients 

of ao and do are constants by taking g[^l = g-^ and r'^l = r— into eq. (fT2|) and eq. p^ . which 

"0 ao 

gives a trivial solution. What is more, we can learn those equation q^ ai + 2 i bo — g^^'x*^! = 0) 
— ai + Vx di — 2 i Co = 0,— g'-'^^d-i + ga_i -|- 60 = and — rt^'o-i + rd-i + cq = from eq. p^ and 
eq. (|14p . Under the condition aidia_id_i 7^ 0,we can get that a_i , ai and di are constants from 
eq. (|12p and eq. (|14p . This means that our assumption does not suppresses the generality of the DT of 
the FL system. Based on eq. (jl2p and (jl4p and above analysis, let Darboux matrix T be the form of 



Here ai,(ii,6o)Co are undetermined function of (x, i), which will be expressed by the eigenfunction 
associated with Ai and A2 in the FL spectral problem, and a^i,d-i are constants. First of all, we 
introduce n eigenfunctions ijjj as 

V'i = ^ ' = l,2,....n,(l)j = (f)j{x,t,Xj), Lfj = (fj{x,t,Xj). (16) 
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Theorem l.The elements of one-fold DT are parameterized by the eigenfunction ipi (i = \,2) asso- 
ciated with Aj, as follows 



Ti(A;Ai;A2) 



aiX + o-iX ^ bo 

Co diX + d-iX^^ 



(17) 



and then the new solutions gt^l and rt^l are given by 
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Proof: Note that a_i and (i_i are constants, which are derived from the eq.(|12p and eq.(|14p. 
respectively. Prom eq. ()14p and transformation eq. ljlSp . new solutions are given by 



a-i , ^0 [11 
d-i a_i 



^cLi ^ Co 
a_i a_i 



(19) 



By using a general fact of the DT, (i.e.), Ti(A; Ai; A2)|A=Ai'0j = 0, i = 1, 2, new solutions are given as 
eq. (jlSp . Further, by using the explicit matrix representation eq. (|17p of Ti, then ■0^^^ is given by 
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After, a tedious calculation shows that Ti in eq. (ll7p and new solutions indeed satisfy eq. (llip and 
eq. (ll3p or (equivalently eq. (ll2p and eq. (ll4p ). So FL spectral problem is covariant under transfor- 
mation Ti in eq. (ll7p . and thus it is the DT of eq.(l2]) and eq.dSj). □ The remaining issue is how to 
guarantee the validity of the reduction condition, i.e., q^^^ = (A^^)* . We shall solve it at the end of this 
section by choosing special eigenfunctions and eigenvalues. 

2.2 N-fold Darboux transformation for FL system 

The key task is to establish the determinant representation of the n-fold DT for FL system in this 
subsection. For this purpose, we set 



D 
A 



a 

d 

b 

c 



a, d are complex functions of x and t 

b, c are complex functions of x and t > , 



as in ref. |53j . 

According to the form of Ti in eq. ijlSp . the n-fold DT should be of the form |53j 

n 

Tn = Tn{X; Xi, X2, • " J ^2n) = PlX\ 



(20) 
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p 



Here P_„ is a constant matrix with a_„ 
Specifically, from algebraic equations, 
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coefficients Pi is solved by Cramer's rule. Thus we get determinant representation of the T„. 
Theorem2. The n-fold DT of the FL system can be expressed by 
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Next, we consider the transformed new solutions (^["^1, rl"^!) of FL system corresponding to the n-fold 
DT. Under covariant requirement of spectral problem of the FL system, the transformed form should 
be 

StV'^i = {jy? + qI'^Ia + vbi"! + y-i^A-^ + ^ ja-2)^n = fNv,n, (24) 

7 



with 



J 



-i 
i 





and then 



-i + pq^' 



Tut + TnV = r„ 



(25) 

Substituting T„ given by eq. (pO|) into eq. (p5|) . and then comparing the coefficients of A^("+^\ it yields 

=(? + 6_(„_i), rW =r + c_(,_i). (26) 



Furthermore, taking c_(„„i) which are obtained from eq. 

are given by 



Jn] 



r + 



(n-l) 



then new solutions (gM^j-N) 

(27) 



with 



(n-l) 



(n-l) 



A2n-l<^2ri-l 
A2n</'2n 

A2n-lV'2n-l 
A2nV'2n 



Ai y^i 
A2 ^2 



A2„„iy'2n-l 
■^2n </'2n 

ArVi 

A2"V2 



A 



2n-l'P2n-l 
n-l. 



("-2) , 
2„-l 'P2n-1 

A2n <P2n 

, -(n-2) 
Ai Vl 
, -(n-2) 
A2 V2 



, -(n-2) 
A2„_i V52n-1 



-Ai "4>1 
-A2"02 

-A 



2„-l</^2n-l 
-A^„"'(;/'2n 



-Ar>i 
-A2 ^2 



^2n 4>2n 



-{n-2) 



'2n 



¥'2n 



"A2n-l'/^2n-l 
^2n 



-^2n^2n 



and the new eigenfunction -f/;^"' of A,fc is 



A^ 




Ar 




■■ A-("-)^. 


^ -(n-l) 

Afc Vfc 




A^ 


^1 


Ar 


-Vl . 


• . A^^-^Vi 


, -(n-l) 
Ai Vl 


Ar>i 


A^^ 


62 


Ar 


'V2 


.. A-("-V2 


, -(n-l) 
A2 V2 


A2>2 



.-{n-2)^ 
^2n-l 



^2n-l 



>^2n^ 



'2n 



\n— 1 
A 



'2n '/^2n 



A 



-(n-2) 
2n S 



^-(n-l) / 
A2n-1 V52n-1 A2„_i02n-1 



02n 



A 



-(n-l) 
2n 



^2n 



^2n^2n 





ArVfe . 


^ -(n-2) 


, -(n- 






A^V^i 


ArVi • 


. -(n-2) 

\^ Vl 


A-("- 


-^Vl 


Ar>i 




ArV2 . 


, -(n-2) 
A2' V2 


> -(n- 
A2 


-^^2 


A2->2 



A2„_i'^2n-1 
A2nV^2n 



V n-l 
'2»i 

A 



A2„-l'?'2n-l 
1 n-l 



2n <?'2n 



, -(n-2) 
A2„_i V'2n-1 
-(n-2). 



^2n-l 



A 



2n 

\Wn 



^2t. 



A 



-{«-!), 
2n 



'2n-l A2^_iV52n-l 



i'2n 



A2n ¥'271 



We are now in a position to consider the reduction of the DT of the FL system so that g'"! = {A""')*, 
then the DT of the FL equation is given. Under the reduction condition q = r*, the eigenfunction 

V'fc = I '^'^ ) associated with eigenvalue Afc has following properties, (fk* = fi, 'Pk* = 4>u ^k* = 
\ V^k J 

where k ^ I. For example, setting I = 1, 3, . . . , 2n — 1, then choosing the n distinct eigenvalues and 
eigenfunctions in n-fold DTs in the following manner: 

A« o V'; = ( ) , andAs; = A;,_i, ^ ^21 = ) (28) 

so that = (r["l)* in eq. ([27j) . Then with these paired-eigenvalues Aj and paired-eigenfunctions 
= 1, 3, . . . , 2n — 1) is reduced to the n-fold DT of the FL equation. Notice that the denominator 
Wn of gl'^l is a modulus of a non-zero complex function under reduction condition, so the new solution 
gil"] is non-singular. 

3. The n-order rogue waves and their determinant forms 

Using the results of DT above, breather solutions of FL equation are generated by assuming a 
periodic seed solution, then we can construct the rogue waves of the FL equation from a Taylor series 
expansion of the breather solutions. 

Set a and c be two complex constants, then q = ceicp iiiax + {- — — ^ is a periodic 

a 

solution of the FL equation, which will be used as a seed solution of the DT. Substituting q = 

cexp(z(ax -|- ( ^'^ — ^ c^)^)) iiito the spectral problem eq.dH) and eq.(l5|), and using the method of 

a 

separation of variables and the superposition principle, the eigenfunction ip2k-i associated with X2k-i 
is given by 

f4>2k-i{x, t, X2k-iJ\fCiT^i{x, t, A2fc-i)[l}+C2ti72(x, t, A2fc-i ) [Ij+Cs^^l (aj, -^2fc-l)[2]+C4TO2(x, t, Xlk_i)[2 
\f2k-i{x, t, \2k-1)) t, A2fc-i)[2}+C2ro2(x, i, A2fc_i)[2}fC3G7^(a;, t, \l^_^)[l}\Ciwl{x, t, ^k-i)^^. 

Here 

2,1 1 r„,-i\2 



.(29) 



■aji{x,t, A2fc-i)[l] 
wi{x,t, A2fc_i)[2] 



W2{x,t, A2fc-l)[l] 
W2{x,t,\2k-l)[2] 



/ ^ 1 2aA2fc_i-F-l , 1 , AoHY \ 

/ exp(Vg(A2fc-i)(- x | \ ty^iiax-a^- '-^^)t)) \ 

a + 2X2k-i^ -WS{X2k-i) , /^TT r.l ,2aA2fc-i^+l ,, I., ^{a¥lf 2^^^ 

I exp V'5(A2fc-i)(-x4— — ^t)--i{axM^ c')t)) j 

\ 2X2k-ica ^ 2 4aA2fc-i 2a/ 

/ exp(-V^(A2fc-i)(-3 ; | '^2 —<?Y)) ^ 

2 4aA2fc_i 2 a 

a + 2A2fc-i' + Vg(A2fc-i) . 1 2oA2fc-i^,. I., ^{ai-lf 2^,^^ 

\ exp A2fc-i) -xH— —t)--i{ax-{{ c^)t)) 

\ 2X2k-ica ^ 2 4aA2fc-i 2a/ 



^ . X \_ I ■CUi{x,t,X2k-l)[l] \ ^ . y ro2(x,t, A2fc_l)[l] 

Wi{X,t,X2k-l) — „ . x ^^r,^ , ^2yx ,t, A2k-l) 



VJi{x,t,X2k~l)[2\ J ' '-^^-''-'-^^-^^ ^ Zn2(x,t,A2fc^l)[2] 

S{X2k-i) = -a^ - 4A2fc-i' + 4A2fc_i2(c2a2 - a)(A; = 1, 2, . . . , n). 

Here a,c,x,t S M, Ci,C2,C^,C4 E C. Note that wi{x,t, X2k-i) and ■uj2{x,t, X2k-i) are two different 
solutions of the spectral problem eq.([4]) and eq.([5]), but we can only get the trivial solutions through 
DT of the FL equation by setting eigenfunction ip2k-i be one of them. What is more, we can get 
richer solutions by using (I29p . 



3.1 The first-order rogue waves generated by first-order breather solutions 
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Under the choice in eq. ()28p with one paired eigenvalue Ai = ai + i(3i and A2 = ai —if3i, the two- fold 
DT eq. (|27p of the FL equation implies a solution 



A2Al(-A292l<^l* + Al0l0l ) 

with 01 and (fi given by eq. (p9|) . For simplicity, under the condition Ci = C2 = C3 = C4 = 1, 



let c = ^ ^""^ so that Im(— — 4A2fc-i'^ + 4A2fc-i^(c^a^ — a)) = 0, we get the first 

a 

order breather q^^^ |46| . Furthermore, by letting a — >• 2(ai^ + /3i^) in (|3U|) with Im(— — 4A2fc_i^ + 
4A2fc-i^(c^a^ — a)) = 0, its first order breather gl^l becomes a rogue wave (7™ [46j . 

3.2 The n-order rogue waves and their determinant forms 

In order to make the higher order rogue waves informative, we modify Ci , C2 , C3 and C4 in the 
equation ([29]) 

fe-i 



Ci = Ko + exp(--i5(A2fc_i) ^ Ji(A2fc-i - Aq)^ 

i=o 

^ fc-i 
C2 = Ko + exp(-iS'(A2fc-i) ^ Jj(A2fc-i - XoY) 

j=0 

^ k-l 

C3 = Ko + exp(--i5(A2fe-i) ^ ii(A2fe-i - Aq)-^ 



2 

fc-i 

C4 = /iTo + exp(-i5(A2fc-i) J2 hi^2k-i - Ao^ ) (31) 

j=0 

Here /Cq) Jj^Lj G C. Note that A2fc-i = Aq = — ac + i—\/—a^c^ + 2a is the zero point of 5'(A2A;-i). 
J J 2 2 

In this way, higher order rogue waves can be constructed from the higher breather solutions. In 
other words, let A2fe-i — s- + + 2a in n-order breather solutions, n-order rogue waves 

can be given. Generally, in comparison to the method of limiting the breather solutions, the method 
of making rational eigenfunction below may be more direct and the rogue wave can be shown in 
determinant forms. 

\ \ 

Substituting eq. (f3T]) into eqs. (f29l) . by assuming A2/C-1 — )• + i-^^ —a^(? + 2a, eigenfunction 

■02^-1 associated with A2A;-i become rational eigenfunction -0^- For simplicity, when a = l,c = —1, 
rational eigenfunction 0^ has the following form 



\{-2 - 2i)Kx - 4Kt + (-1 + i){-Ll + + 2i + 2iKo - 2iLo + 2Lq)) exp(^i(x + 3t)) ' 

1 

((2 + 2i)Kx + AKt + i){-2i - + - 2iKo - 2iJo + 2Jo)) exp(--i(x + 3t)) 

K = -Lo + l-i-iKo + Ko + Jo (32) 
Substituting eigenfunctions eq. (p2]) into eqs. (p7j) . we can get the first order rogue wave q!nl in the 



form of determinant. The dynamical evolution of q^rw for the parametric choice Kq = 1, Jq = Lq = 10 
is plotted in the Figure 1, which control the position of the first-order rogue wave by choosing the 
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parameters Kq.Jq and La- Similarly, the corresponding density plot is portrayed in the Figure 2. 
Theorem 3. For the n-fold DT, the n-order rogue wave qlw is of the form 

q + 



drw 



-(n-1) 



(33) 



with 



"'n2 



un * 
"■n2 



'''n-12 



V-ll 



''■n-22 



'^n-2l 



'''n-32 



1^. 



-(n-1) 



"-712 



* 

"'n2 



"'n-12 
'''n-11 



'n-22 
,1 

■-71-21 



'V-32 



hi 



'■n-12 



''n-11 



n-22 



''ri-21 



''n-22 



n-31 



'''n-32 
Ll * 
''n-31 



'n-32 



''n-31 



The final form of Tj.„(A) has the form, 



Trn — Trn{^', Ai, A2, ■ ■ ■ , A2n) 



ui 



(n-2)l 

* 

-(n-2)2 



(n-l)2 

* 

(n-l)l 



-(n-2)l 
I * 
-(n-2)2 



(n-2)l 
(n-2)2 

I 

-(n-2)l 
(n-2)2 



un 

"'-(n-l)2 
'^-{n-l)l 



nl 



-h 



-n2 



-nl 



-n2 



mill 



rn 12 



Wrn 
{Trn) 21 



Wrn 
{Trn) 22 



\ Wr. 



Wr. 



(34) 



with 



(^rn) 11 



{T. 



rn) 12 



{Trn) 



21 



'''nl 

ul * 
"'n2 



un 
"nl 
un * 
"■n2 



"-nl 
ul * 
"'n2 



"-nl 
un * 
"■n2 



^n2 
ul * 
"nl 



"'n2 
un * 
"nl 



'^-12 

ul * 
"n-11 



un 

'''n-12 
un * 
"n-11 

Xn-1 

"n-12 
ul 

"n-11 



h: 
K: 



n-12 



n-11 

y^n-l 

"n-11 
ul 

'''n-12 



'n-ll 
n * 
n-12 



A— 2 

"n-21 
ul * 
"'n-22 



un 

'''n-21 
un * 
'''n-22 



"n-21 
ul 

"'n-22 



"'n-21 
un * 
"'n-22 



'''n-22 
ul 

"'n-21 



"'n-22 
un * 
"n-21 



'''n-32 
ul * 
"n-31 



"'n-32 
un * 
'''n-31 

Xn-3 

'''n-32 
ul 

"'n-31 



"'n-32 
un * 
'''n-31 

Xn-3 

"n-31 
ul 

"n-32 



in 

"■n-31 
un * 
'''n-32 
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X-in-2) 



-(n-2)l 
* 

(n-2)2 



(n-2)l 
1-: 





un 

''■-(n-2)2 



"'-(n-2)l 
ul 

"'-{n-2)2 



-(n-2)l 
I * 
-(n-2)2 





(n-2)2 
(n-2)l 



-(n-2)2 
I * 
-(n-2)l 





,1 

''-(n-l)2 
1 * 
-(n-l)l 



"-(n-l)2 
un * 
'^-(n-l)l 

"-(n-l)2 
ul 

'^-{n-l)l 



h 
h 



-(n-l)2 

(n-l)l 

A-("-i) 



-(n-l)l 
(n-l)2 



-(n-l)l 
I * 
-(n-l)2 



A-" 

"-nl 
ul * 
"-n2 



"--nl 
un * 
"-n2 



nl 



n2 



''-nl 



-n2 



n2 



nl 



''-n2 



-nl 



Here 



rn 22 



hi2 

ul * 
"'nl 



''■n2 
i,n * 
''■nl 



''-ri-ll 
ul * 
"n-12 



/ 



un 
""n- 



12 



yn-2 

''"n-22 
ul * 
"n-21 



ul * 
"'n-32 



"'n-22 "-n-Sl 



''"n-21 '''n-32 



X-(n-2) 







' -(n-2)2 "--(n-l)! 
' -(n-2)l "-(n-l)2 



"'-(n-2)2 '^-(n-l)l 
'''-(n-2)l '''-(n-l)2 



^1 



n2 



nl 



''-n.2 



-nl 



d 

4i = ^((Ao + <J)"'0i(Ai = Ao + (5))|5=o,m 



-n, 



[n 



I),--- ,-l,0,l,--- ,n 



1, n, Z 



1,2, 



2n. 



hL2 = ^((^0 + 5)"'</'i(Ai = Ao + (5))|5=o,m = -n, -(n - 1), • • • ,-1,0,1,- •• , n - 1, n, / = 1, 2, • • • 2n. 

Case 1). When n = 2, substituting eq. (j29|) into eq. ([33|) . one can construct the second-order rogue 
waves with seven free parameters. Note that under the condition Ji >> Jq and Li >> Lq except 
for Jo = —Lq,Ji = — Li, the second-rogue can spht into three first-order rogue wave (triplets rogue 
wave) [41] rather than two. The dynamical evolution of Igi^p for the parametric choice Kq = 1, Jq = 
— Lq, Ji = Li,Lq = — 5,Li = 20,a = l,c = — 1 is plotted in the Figure. 3 and the corresponding 
density plot is shown in the Figure. 4. There is another kind second-order rogue wave, for example, 
is higher than second-rogue above, which is a a fundamental pattern [22]. The dynamical 
evolution of If/i^P for the parametric choice Jq = —Lq,Ji = —Li,a = l,c = — 1 is plotted in the 

Figure. 5 and the corresponding density map is portrayed in Figure. 6. Note that has only two 

free parameters a and c under the condition Jq = —Lq, Ji = —Li. 

Case 2). When n = 3, substituting eq. (f29]) into eq. (f33]) can lead to the third-order rogue waves with 
nine free parameters. Note that under the condition J2 >> Ji, L2 >> Li{i = 0, 1) or Ji >> Jj, Li » 
Li{i = 0,2) except for Jq = —Lq,Ji = — Li, J2 = —L2, the third-rogue can split into six first-order 
rogue wave rather. Circular rogue wave [S] may be constructed by the condition J2 >> Ji,L2 >> 
Li{i = 0,1) except for Jq = —Lq,Ji = — Li, J2 = —L2. The dynamical evolution of | g'r'E P (circular 
pattern) for the parametric choice Kq = 1, Jq = —Lq,Ji = — Li, J2 = -^2,-^0 = 5, Li = 5,L2 = 
3000, a = l,c = — 1 is plotted in the Figure. 7 and the corresponding density plot in Figure. 8. At 
the same time, triplets rogue wave may be constructed by the condition Ji >> Jj, Li » L-iii = 0, 2) 
except for Jq = —Lq, Ji = — Li, J2 = —L2. The dynamical evolution of | q'I.'^ p (triangular pattern) for 
the parametric choice Kq = 1, Jq = —Lq, Ji = Li, J2 = —L2,Lq = 0, Li = 500, L2 = 0, a = 1, c = — 1 
is plotted in the Figure. 9 and its density map in Figure. 10. Similarly, there is another kind of 
third-order rogue wave, for example, is higher than third-rogue above. The dynamical evolution 

of I p (fundament pattern) for the parametric choice Jq = —Lq, Ji = —Li, J2 = —L2, a = 1, c = — 1 
is plotted in the Figure. 11. Similarly, the density plot of Figure. 11 and the corresponding density 
plot in Figure. 12. Note that has only two free parameters a and c under the condition 

Jo = —Lq, Ji = — Li, J2 = —L2. 

Case 3). When n = 4, substituting eq. ([29]) into eq. (p3]) one can generate the four-order rogue 
waves with eleven free parameters. Note that under the condition J3 >> Ji,L3 >> Lj(i = 0,1,2) 
or J2 >> Ji,L2 >> Li{i = 0,1,3) or Ji >> Ji,Li » Li{i = 0,2,3) except for Jq = — Lq, Ji = 
—Li, J2 = —L2, J3 = —L3, the four-rogue waves can split into ten first-order rogue wave rather. One 
kind of circular rogue wave [H] may be constructed by the condition J3 >> Jj,L3 >> Li{i = 0, 1,2) 

except for Jq = — Lq, Ji = — Li, J2 = — L2, J3 = —L3. The dynamical evolution of Igitlp, which 
is a fourth order rouge wave consisting of a ring structure (outer, seven peaks) and a fundamental 
pattern of the second rogue wave (inner), for the parametric choice Kq = 1, Jq = —Lq, Ji = — Li, J2 = 
— L2, J3 = L3,Lq = Li = L2 = 5,1/3 = 2000, a = l,c = — 1 is plotted in the Figure. 13 and the 

12 



corresponding density plot in Figure. 14. In addition, the dynamical evolution of |(7r«)P, which is 
another circular pattern rogue wave [44], for the parametric choice Kq = 1, Jq = Lq, Ji = —Li, J2 = 
—L2,J3 = L^jLq = 50, Li = 10, L2 = 0, L3 = 2000, a = l,c = — 1 is plotted in the Figure. 15. 
and its density map is portrayed in Figure. 16. Note that the inner structure of Figure 15 (or 16) 
is a triangular pattern of a second order rogue wave. At the same time, the triangular pattern of 
the rogue wave may be constructed by the condition Ji >> Ji,Li » Li{i = 0,2,3) except for 
Jo = —Lq,Ji = — Li, J2 = —L2,J3 = —L3. The dynamical evolution of | gitl P (triangular pattern) 
for the parametric choice Kq = 1, Jq = — Lq, Ji = ^1,^2 = —L2,J3 = —L^^Lq = 0, Li = 2000, L2 = 
1/3 = 0, a = 1, c = — 1 is plotted in the Figure. 17 and its density plot in Figure. 18. A pentagon 
pattern of the rogue wave may be constructed by the condition J2 » Ji,L2 » Li{i = 0,1,3) 

except for Jq = — Lq, Ji = — Li, J2 = — -L2, J3 = —L3. The dynamical evolution of jg'itlp (pentagon 
pattern) for the parametric choice Kq = 1, Jq = — Lq, Ji = — Li, J2 = L2, J3 = —L^^Lq = Li = 
0, L2 = 2000, L3 = 0, a = 1, c = — 1 is plotted in the Figure. 19. Similarly, the density plot of Figure. 
19 is correspondingly shown in Figure. 20. Similarly, there is another kind four-order rogue wave, 
for example, |(7rtlp is higher than four-rogue above. The dynamical evolution of Igl^p (fundamental 
pattern) for the parametric choice Jq = —Lq, Ji = — Li, J2 = —L2, J3 = —L3, a = 1, c = — 1 is plotted 
in the Figure. 21 and the corresponding density plot is portrayed in Figure. 22. Note that l^rtlp 
only has two free parameters a and c under the condition Jq = —Lq, Ji = — Li, J2 = —L2, J3 = —L3. 
According to the above analysis, the n-order rogue waves may be controlled by 2n + 3 free parameters. 

4. Conclusions 

Thus, in this paper, considering FL system of equation which describes nonlinear pulse propagation 
through single mode optical fiber, the determinant representation of the N-fold DT for the FL system 
is given in eqs. (i23]) . The n-order rogue wave in eq. (f33]) of the FL equation is obtained by this repre- 
sentation from a periodic seed under the reduction condition eq. (l28p . Several interesting patterns of 
the rogue wave are plotted by choosing suitable parameters. 

Our results provide an alternative possibility to observe rogue waves in optical system. Moreover, 
from the one- fold DT, it is interesting to observe that the DT of the FL system exhibits the following 
novelty in comparison with other integrable models like the AKNS and the KN systems: the DT 
matrix of the FL system has three different terms depending on A. Thus, the DT as well as the rogue 
wave of the FL system holds novel features, in comparison with the DT and rogue wave solutions of 
the standard integrable counterparts like the AKNS and the KN systems. The n-order rogue waves of 
the FL equation are constructed by using the determinant representation of the DT with 2n + 3 free 
parameters. 
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Figure 1. The dynamical evo- 
lution of IgiwP with specific pa- 
rameters Kq = 1, Jo = Lq = 
10, a = l,c = -1. 



Figure 2. Contour plot of the 
wave amplitudes of [gi^p for the 
values used in Figure 1. 
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Figure 3. The dynamical evo- 
lution of with specific 
parameters Kq = 1, Jq = 
—Lo,Ji = Li,Lq = — 5,Li = 
20, a = l,c = -1. 



Figure 4. Contour plot of the 
wave amplitudes of \q?w\'^ for the 
values used in Figure 3. 




Figure 5. The dynamical evo- 
lution of with specific 
parameters Jq = —Lq, Ji = 
—Li, a = 1, c = —1. 




Figure 7. The dynamical evo- 
lution of with specific 
parameters Kq = 1 , Jq = 
—Lo,Ji = —Li,J2 = L2,Lq = 
5,Li = 5,L2 = 3000, a = l,c = 
-1. 



Figure 6. Contour plot of the 
wave amplitudes of jgi^p for the 
values used in Figure 5. 
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Figure 8. Contour plot of the 

wave amplitudes of Igriip for the 
values used in Figure 7. 
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Figure 9. The dynamical evo- 
lution of with specific 
parameters Kq = 1 , Jq = 
—Lo,Ji = Li,J2 = —L2,Lq = 
0,Li = 500, L2 = 0,a = l,c = 
-1. 




Figure 1 1 . The dynamical evo- 
lution of with specific 
parameters Jq = —Lq, Ji = 
—Li, J2 = —L2, a = 1, c = — 1. 




Figure 13. The dynamical evo- 
lution of Igitlp with specific 
parameters Kq = 1 , Jo = 
—Lq,Ji = — Li, J2 = — ^2,^3 = 
LsjLq = Li = L2 = 5,L3 = 
2000, a = l,c = -1. 



Figure 10. Contour plot of the 

wave amplitudes of | qlw \ ^ for the 
values used in Figure 9. 




Figure 12. Contour plot of the 
wave amplitudes of for the 

values used in Figure 11. 




Figure 14. Contour plot of the 
wave amplitudes of Igitlp for the 
values used in Figure 13. 
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Figure 15. The dynamical evo- 
lution of l^rap with specific pa- 
rameters Kq = 1, Jo = Lq,Ji = 

— Li,J2 = — L2,J3 = L^jLq = 

50, Li = 10, La = 0,^3 = 
2000, a = l,c = -1. 
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Figure 16. Contour plot of the 
wave amplitudes of for the 

values used in Figure 15. 




Figure 17. The dynamical evo- 
lution of Igitlp with specific 
parameters Kq = 1, Jq = 
—Lq,Ji = Li,J2 = —L2,J3 = 
-L3,Lo = 0,Li = 2000, La = 
L3 = 0, a = l,c = — 1. 
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Figure 18. Contour plot of the 
wave amplitudes of Igitlp for the 
values used in Figure 17. 
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Figure 19. The dynamical evo- 
lution of with specific 
parameters Kq = 1 , Jo = 
—Lq,Ji = — Li,J2 = ^2,^3 = 
~Ls,Lq = Li = 0,^2 = 
2000, L3 = 0,a = l,c = -1. 
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Figure 2 1 . The dynamical evo- 
lution of Igitlp with specific 
parameters Jq = —Lq, Ji = 
—Li,J2 = —L2,J3 = —L3,a = 
l,c= -1. 



6 
4 

2 

t a 

-2 
-* 
s 



Figure 20. Contour plot of the 

wave amplitudes of for the 

values used in Figure 19. 
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Figure 22. Contour plot of the 
wave amplitudes of for the 

values used in Figure 21. 



19 



